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1. Introduction 

Lattice QCD simulations are currently being performed with two dynamical degenerate light 
quarks with a mass close to their physical value as well as the strange quark using a number of 
different discretization schemes with the most common being Wilson-improved, staggered and 
chiral fermions. Furthermore, simulations at several lattice spacings and volumes are becoming 
available enabling a comprehensive study of lattice artifacts. 

The focus of this contribution is the evaluation of form factors (FFs) and moments of parton 
distributions of the nucleon, which are being measured in many experiments. The characterization 
of nucleon structure is considered a milestone in hadronic physics and experiments on nucleon FFs 
started in the 50s. A new generation of experiments using polarized beams and targets are yielding 
high precision data spanning a larger range of momentum transfers. FFs provide ideal probes of 
the charge and magnetization densities of the hadron as well as a determination of its shape. 

Lattice techniques to extract nucleon matrix elements connected to FFs and moments of gener- 
alized parton distributions (GPDs) are well developed. The connected diagram where on operator 
couples to a valence quark is straightforward to compute and has been evaluated by a number of 
lattice groups [1]. For iso-vector nucleon matrix elements, in the isospin limit, this is the only 
contribution. Like most collaborations, we use non-perturbative renormalization of these matrix 
elements and furthermore we subtract €?(a 2 ) -terms computed perturbatively to improve the ex- 
traction of the renormalization constants [2]. Using simulations of Nf = 2 twisted mass fermions 
(TMF) at three values of the lattice spacing and different volumes we study cut-off and volume 
effects. We use the lattice spacing determined from the nucleon mass to convert to physical units. 
Heavy baryon chiral perturbation theory (HB^PT) is used to extrapolate lattice results obtained for 
pion masses in the range of about 260 MeV to 470 MeV to the physical point. 

2. Nucleon form factors 

The nucleon matrix element of the electromagnetic (EM) current, = \j/(x)y^Y( x )^ * s written 
in the form Un{p' ,s') y^F\ (q 2 ) + ^(g 2 ) uj^{p,s), where the Dirac F\ and Pauli F2 FFs are 

related to the electric and magnetic Sachs FFs with the relations: Gpiq 2 ) = Fi(q 2 ) — p&i^K? ) 
and GM(q 2 ) = F\{q 2 ) +F2(q 2 ). For the axial vector current, = y/(;e)y^y5yi//(x), the nucleon 
matrix element is of the form un(p', s') y^/sGU^ 2 ) + ^^Gp^ 2 ) \un(p,s). 

The axial charge is well known experimentally. Since it is determined at Q 2 = there is no 
ambiguity associated with fitting the <2 2 -dependence of the FF. As can be seen in Fig. 1, where we 
show recent lattice results using TMF, domain wall fermions (DWF) and a hybrid action of DWF 
valence on staggered sea quarks, there is a nice agreement among different lattice discretizations 
and no significant dependence on the quark mass down to about m n = 270 MeV. 

We take the continuum limit of TMF results by InAtting to a constant, after checking that a 
lineal - inA_t at two values of the pion mass yields consistent results with the constant inA_t. Volume 
corrections are taken into account following Ref. [3]. The volume corrected continuum results are 
shown in Fig. 1. Chiral extrapolation using HB^PT with three fit parameters in the small scale 
expansion (SSE) [4] produces a value of gA = 1.12(8) at the physical point, which is lower than 
the experimental value by about a standard deviation. The large error band is due to the large 
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Figure 1: Left: Lattice data on g A using N F =2 TMF [1], N F =2 + 1 DWF [5] and N F = 2+1 using 
DWF on staggered sea quarks [6]. The physical point is shown by the asterisk. Right: Volume corrected 
continuum TMF results together with the band obtained using HB^PT. 
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Figure 2: Left: Isovector electric and magnetic nucleon FFs at m n ~ 300 MeV using TMF [1], DWF [7], 
hybrid [6] and Clover [8]. Experimental data are shown with the filled circles accompanied with Kelly's 
parametrization shown with the dashed line. Right: Axial nucleon FFs. The solid line is a dipole fit to 
experimental data for Ga(Q 2 ) combined with pion pole dominance to get the solid curve shown for G P (Q 2 ). 



correlations between the A axial charge #aa and the counter-term involved in the chiral expansion. 
A lattice determination of gAA [9] will therefore allow a more controlled chiral extrapolation. 

In Fig. 2 we show recent lattice results on the EM isovector FFs Ge{Q 2 ) and Gm{Q 2 ) at m n ~ 
300 MeV, where we see a nice agreement for Ge{Q 2 ) but a clear disagreement with experiment, 
with lattice data showing a weaker <2 2 -dependence. For Gm{Q 2 ) there is some spread in the results 
that needs to be investigated. For the axial Ga(Q 2 ) FF there is good agreement among TMF, DWF 
and hybrid results. Like in the case of the EM form factors, the <2 2 -dependence of Ga{Q 2 ) is milder 
than what is observed experimentally. TMF results on G P (Q 2 ) at m % ~ 300 MeV on a lattice with 
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spatial L = 2.8 fm and results using a hybrid action of DWF on staggered sea at m n ~ 350 MeV and 
L = 3.5 fm show discrepancies at low Q 2 , that may be due to the smaller volume in our calculations. 

3. Nucleon moments 

In Figs. 3 we compare recent results from ETMC [1], RBC-UKQCD [10], QCDSF [11] and 
LHPC [6] on the spin-independent and helicity quark distributions. All collaborations except LHPC 
use non-perturbative renormalization constants. The ETMC has, in addition, subtracted &(a ) 
terms perturbatively [2]. There is a spread among lattice results with results obtained with the 
hybrid action being lower than those from ETMC and QCDSF. Using HB^PT [12] we extrapolate 
results on A 20 and A20 to the physical point as shown by the curves in Fig. 3. Our estimates for both 
^20 = (*)«-rf an d A20 = (x) Au- Ad are considerably higher than the experimental values. The very 
recent result by QCDSF [1 1] at m % ~ 170 MeV remains higher than experiment and highlights the 
need to understand such deviations. 
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Figure 3: Recent results on {x) u -d and (x) A „_ A ^. The fit is done using HBjPT [12]. 
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